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Cosmic ray air showers are known to emit pulsed radio emission which can be understood as coherent geosyn- 
chrotron radiation arising from the deflection of electron-positron pairs in the earth's magnetic field. Here, 
we present simulations carried out with an improved version of our Monte Carlo code for the calculation of 
geosynchrotron radiation. Replacing the formerly analytically parametrised longitudinal air shower develop- 
ment with CORSIKA-generated longitudinal profiles, we study the radio flux variations arising from inherent 
fluctuations between individual air showers. Additionally, we quantify the dependence of the radio emission 
on the nature of the primary particle by comparing the emission generated by proton- and iron-induced show- 
ers. This is only the first step in the incorporation of a more realistic air shower model into our Monte Carlo 
code. The inclusion of highly realistic CORSIKA-based particle energy, momentum and spatial distributions 
together with an analytical treatment of ionisation losses will soon allow simulations of the radio emission 
with unprecedented detail. This is especially important to assess the emission strengths at large radial dis- 
tances, needed for event-to-event comparisons of the radio signals measured by LOPES in conjunction with 
KASCADE-Grande and for considerations regarding large arrays of radio antennas intended to measure the 
radio emission from ultra-high energy cosmic rays, as with LOFAR or in the Pierre Auger Observatory. 
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1. Introduction 

As part of LOPES (Q, a digital radio interferometer working in conjunction with the KASCADE-Grande ex- 
periment, we have developed a model of the emission mechanism responsible for the radio pulses produced 
by extensive air showers. Describing the emission as coherent geosynchrotron radiation arising from the de- 
flection of electron-positron pairs in the earth's magnetic field, we have performed analytical calculations |2| 
followed by detailed Monte Carlo simulations [3|, which were then used to derive important dependences of 
the radio emission on specific air shower parameters |4|. For the first time, our simulations have employed a 
realistic model of the air shower producing the radio emission. Important air shower properties such as the lon- 
gitudinal ("arrival time") and lateral particle distributions in the shower pancake, the particle energy and track 
length distributions and the overall longitudinal evolution of the shower were taken into account via "standard" 
analytical parametrisations such as NKG functions (lateral distributions) and Greisen parametrisations (longi- 
tudinal evolution). While a number of important results have already been obtained with this model |4), the 
next step in improving the simulations is to substitute the underlying air shower model by a more realistic one 
based on histogrammed distributions derived from CORSIKA [5 1 simulations. 

Here, we present a first result derived with these improved simulations, namely the influence of the air showers' 
realistically simulated longitudinal profiles. These realistic CORSIKA-generated profiles allow us to evaluate 
shower-to-shower fluctuations in the radio emission as well as the influence of the nature of the primary particle 
on the radio signal. 



2. Discussion 

To assess the relative fluctuations in the radio signal between individual air showers, we have simulated 10 
very weakly thinned proton induced air showers and 10 very weakly thinned iron induced air showers with a 
primary particle energy of 10 16 eV and vertical incidence. The geomagnetic field was adopted the same as in 
El . namely 70° inclined with a strength of 0.5 Gauss, which approximately corresponds to the geomagnetic 
field configuration in Central Europe. CORSIKA was used with the GHEISHA2002 model |6| for low-energy 
hadronic interactions and the QGSJET01 model |7| for high-energy hadronic interactions. The relatively low 
primary particle energy was chosen to minimise the computation time in this early stage of analysis. We will 
reconsider the effects for higher energies at a later time in conjunction with a much broader analysis of the 
improvements achieved with the adoption of the other CORSIKA-based particle distributions. 

For the proton and iron showers each, the 10 simulations were averaged to obtain an "average air shower" with 
which one can then compare the individual showers. For this purpose, the individual showers were reviewed by 
eye and the cases with longitudinal profiles differing most from the average profile were selected. Simulations 
of the associated radio signal were then performed with our Monte Carlo code. Only the longitudinal profile 
of the shower development was taken from the CORSIKA simulations. All other air shower properties (energy 
distributions, spatial particle distributions, ...) were retained using the analytical distributions employed in |4]. 
This procedure allows a simple first estimate of the importance of the fluctuations in the shower development. 
For comparison, we also calculated the radio emission from a fully parametrised shower, retaining a Greisen 
parametrisation for the longitudinal development. 

The left panel of figure^shows a comparison of the longitudinal profiles taken into account for the analysis of 
the proton induced showers. In particular, the analytical Greisen parametrisation used so far in the code is very 
similar to the profile of the averaged proton shower. It does, however, overestimate the number of particles 
generated by ~ 15%, as it describes a purely electromagnetic shower of the given energy without significant 
"losses" to a muonic component. The three selected individual shower profiles represent cases of a shower 
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Figure 1. Left: Longitudinal development profiles of vertical CORSIKA-simulated 10 eV proton-induced air showers in 
comparison with an averaged profile derived from CORSIKA and a Greisen-parametrised profile. Right: Simulated radial 
dependence of the 10 MHz radio signal as calculated using the longitudinal profiles displayed on the left. 

developing to its maximum far away from the ground, a shower developing to its maximum close to the ground, 
and a shower having a generally much broader maximum. The right panel of figure^shows the 10 MHz radio 
signal obtained from these longitudinal shower profiles. Due to the overestimation of the number of particles 
in the Greisen parametrisation, the emission is about ~ 15% higher for the fully parametrised shower. Apart 
from this constant enhancement factor, the radial profile of the radio emission is very similar to that of the 
averaged proton shower, which indeed has a very similar longitudinal evolution. The individual showers, in 
comparison, show significant deviations from the radial profile generated by the average shower. While the 
overall emission level (directly related to the number of electrons and positrons) stays similar, the slope of the 
radial emission profile changes. The shower with the maximum close to the ground has a steeper slope, the 
one with the maximum far away from the ground has a flatter slope. This dependence of the radio emission 
on the depth of the shower maximum was already shown in |4|. The shower with the broader maximum has 
slightly lower radio emission over the radial distances considered here. This arises from the lower total number 
of electrons and positrons sampled in the atmosphere's ~ 1000 g cm -2 over the ground, after which the profile 
is "truncated". Overall, the deviations in the field strengths are of the order ~ 30%. 

For the iron-induced showers, the fluctuations in the longitudinal shower development are much smaller, as 
expected from the simplified "superposition model" interpreting an iron-induced shower as a superposition 
(and thus averaging) of 56 nucleon-induced showers. Correspondingly, the deviations of the radio signal 
between the "average shower" and the individual showers are much smaller. This is shown in figure [2] in 
comparison with the fluctuations for the proton-induced showers. As the iron-induced showers grow to their 
maxima at larger distances from the observer than the proton induced showers, the radial profiles of their radio 
emission are somewhat flatter. Additionally, the emission amplitudes are slightly lower. This is due to the fact 
that for iron-induced showers, the fraction of energy going into the muonic shower component is higher (and 
thus the fraction going into the electromagnetic component lower) than in proton-induced showers. 

3. Conclusions 

Using CORSIKA-generated air shower development profiles, we have analysed the influence of shower-to- 
shower fluctuations on the associated radio signal. While vertical 10 16 eV showers induced by proton primaries 
show fluctuations of the order 30%, the fluctuations for iron-induced showers with the same parameters are 
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Figure 2. Radial dependence of the 10 MHz radio emission from vertical 10 eV air showers induced by protons and 
iron nuclei, respectively, comparing the emission from an averaged shower with that from individual showers. 

negligible. Two characteristics of the shower profiles have a direct impact on the radio emission: the distance of 
the shower maximum from the observer on the ground, directly influencing the slope of the radial dependence, 
and the total number of electrons and positrons integrated over the shower evolution, directly determining the 
absolute strength of the emission. Differences in the radio emission of individual showers and also differences 
in the emission from showers induced by different primaries can be related directly to these two effects. 

The analysis we performed constitutes a first step towards composition studies of cosmic rays with radio tech- 
niques. However, it was carried out with only limited statistics at relatively low energies, and only for vertical 
air showers. A more complete analysis at higher energies and various zenith angles will be performed soon in 
conjunction with a detailed analysis of the effects introduced by substituting other analytical parametrisations 
such as those describing energy and pitch-angle distributions by more realistic CORSIKA-derived data. It is, 
however, clear that the influence of differences in the longitudinal shower profiles on the radio emission will 
still mainly be related to the spatial distance of the shower maximum to the ground and the total number of 
electrons and positrons integrated over the shower evolution. 
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